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Abstract. In ad hoc wireless networking, units are connected to each other rather than to a central, fixed, 
infrastructure. Constructing and maintaining such networks create several trade-off problems between 
robustness, communication speed, power consumption, etc., that bridges engineering, computer science 
and the physics of complex systems. In this work, we address the role of mobility patterns of the agents 
on the optimal tuning of a small-world type network construction method. By this method, the network 
is updated periodically and held static between the updates. We investigate the optimal updating times 
for different scenarios of the movement of agents (modeling, for example, the fat-tailed trip distances, 
and periodicities, of human travel). We find that these mobility patterns affect the power consumption in 
non-trivial ways and discuss how these effects can best be handled. 

PACS. 89.75.He Networks and genealogical trees - 89.70.-a Information and communication theory 


1 Introduction 

We are living in an era with an increasing amount of wire¬ 
lessly linked devices. These could be mobile phones (or 
other devices for human communication), environmental 
sensors, military appliances, etc. There is an exciting pos¬ 
sibility of linking these in a decentralized, self-organized 
way, without a central control or a pre-assigned connec¬ 
tion topology and without the need of building a central¬ 
ized infrastructure (which takes time, money and plan¬ 
ning). Such mobile ad hoc networks (often acronymed 
“MANETs”), inspired by distributed systems in nature, 
could potentially be more robust and energy efficient than 
traditional architectures m- Needless to say, in disasters, 
military conflicts and other rapidly changing conditions, 
they would be very useful. Despite this promising outlook, 
ad hoc networks are not yet in widespread use—for one 
reason because it is hard to engineer a system to be self or¬ 
ganized. The challenge to design ad hoc networks is an in¬ 
terdisciplinary problem uniting engineering and computer 
science with the natural sciences In this work, 

we will use a physics approach to study a problem in the 
interface of ad hoc networks and human mobility. 

An ad hoc network needs to simultaneously be robust 
to failures, efficient in transferring information and con¬ 
sume as little energy as possible. The energy consumption 
is controlled by letting only node pairs conneted by a link 
relay packets between them. The network topology should 
also be fair to the agents, so not some agents need to con¬ 
sume much more energy than others. To reduce the energy 
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consumption, the network does not have to relay packets 
from all other devices within range, a reduced number of 
them could be enough mm- The power consumption P 
of a radio link increases rather fast with the distance r be¬ 
tween the devices. In the literature, this relationship has 
been modeled like with 6 ranging between 2 and 4 [T2] . 
Note that there is also a cost associated with receiving a 
packet (that is independent of whether, or not, there is a 
link to the sender), that we, and most of the literature [15] , 
ignores. The motivation is that this cost is lower and more 
uniformly distributed than the sending cost. 

Given its energy budget, an agent in an ad hoc net¬ 
work would then have a choice to have many short-range 
connections, or a few long-range connections. More short- 
range connections would make the network denser and 
more robust against failures, but it would also decrease 
the communication efficiency and reliability of the net¬ 
work. The reason is that an information packet then needs 
to travel more steps to reach from a source to a target. 
The needed energy to route the packet and the chance of 
passing a failing device will increase with increasing link 
density. The good news is that trade-offs need not cost 
that much. This insight comes from the highly cited work 
by Watts and Strogatz on the small-world network phe¬ 
nomenon m- In their model, they start from a locally 
ordered network and show that you only need a small 
fraction of long-range connections to obtain a large de¬ 
crease in the average distance (number of steps from ran¬ 
dom pairs of source and target). Inspired by this. Holme, 
Kim and Fodor (HKF) proposed a method to construct 
the connection topology of an ad hoc network |7] . In their 
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method, each agent has a window of power consumption 
[-Pmim -Pmax] and they commit to the network to spend 
power in that range. Then they fill up their short-range 
connections except occasionally (controlled by a param¬ 
eter) when they connect to a more distant agent. (HKF 
also consider a version of this model where the agents are 
separated into those only making short- and those only 
making long-range connections, with similar properties. 
In this work, we only study the former version.) However, 
HKF assume that agents are static in space. Even though 
there are such applications (e.g. environmental sensors), 
this assumption limits the application rather severely. 

The HKF protocol can quite straightforwardly be ex¬ 
tended to handle mobile agents. One could simply leave 
the topology intact for some period of time and then re¬ 
configure it according to the protocol. However, in that 
case, one faces another trade-off problem. Each time one 
reorganizes the network one needs to spend energy, both 
for numerical processing and communication. So reorga¬ 
nizing too often is not preferable. Reorganizing too rarely 
is also not an option since the movement of the agents 
will deteriorate the original optimization—the energy con¬ 
sumption will increase and there will be a heightened 
chance of links breaking because devices go out of range 
from each other. In this paper, we investigate this trade-off 
problem—how often should one reconfigure the network 
and how does this depend on the structure of the move¬ 
ment patterns. There are many models for human mobility 
patterns, ranging from parsimonious to elaborate [2]. We 
take the former approach and investigate three stylized 
models capturing some simple empirically verified mobil¬ 
ity patterns. One uniformly random scenario, mostly as 
a reference model; one Levy flight scenario (where agents 
take random steps of power-law distributed length); and 
one pattern mimicking periodic commuters. Both the last 
two scenarios are claimed to be typical in human mobil¬ 
ity [HISEE]. 

In the remainder of this paper, we will go through the 
precise definition of the protocol and the mobility models; 
present the numerical results and, finally, put these in a 
wider perspective. 

2 Simulation preliminaries 

2.1 Overview of the simulation 

Our problem is, as described in the Introduction, by na¬ 
ture very complex. Still we try to keep the simulation as 
simple as possible, at the expense of some realism. That 
much said, we do believe some of our insights extends to 
real routing protocols. Our first simplifying assumption 
is that the network is maintained no matter where the 
agents are located—if they are far away, the connection 
just costs more power. As we do not want to relax too 
many assumptions of Ref. [7], we let our model inherit 
this aspect. In a further study, it would be natural to im¬ 
pose a restricted connection range. We also assume agents 
never join or leave the system. Probably the conclusions 
would be qualitatively the same of we had a turnover of 


agents (usually called “churn” in the engineering litera¬ 
ture). We also assume the agents are confined to the unit 
square. It is reasonable to assume a practical ad hoc net¬ 
work is confined to a limited region of space |I0j , then the 
particular shape of the region should not matter much as 
long as the connection range of the devices is longer than 
the longest distance within it. 

Furthermore, we assume the network of radio links are 
updated at regular intervals (with a time Du between each 
update). This set up is similar to real implementations of 
ad hoc networks m- At each update, the connection net¬ 
work is re-optimized according to the HKF protocol [7] 
(discussed in detail below) and then kept constant until 
the next update. Throughout the simulations the agents 
are moving according to a mobility model. We discuss 
these models below. 

2.2 The Holme-Kim-Fodor model 

Given a static point pattern of agents, and the further 
assumptions mentioned above, we construct a network of 
radio connection by the Holme-Kim-Fodor model [7]. As¬ 
sume an agent i has position {xi,yi) and is connected to 
ki other agents Tj. Then the power consumption of i is 

~ ’ ( 1 ) 

where the exponent S captures the distance dependence 
of the radio device. 5-values in the range [2,4] has been 
used in the literature. Following Ref. [7], we use 6 = 2. 
This equation neglects the cost of processing incoming 
packets. This is standard in the literature m and can 
be motivated by it contributing less to the energy cost 
than sending packets, and also more equal for different 
agents (especially in a setting with mobile agents). 

For fairness, we assume an interval [Pinini T^max] of power 
consumption such that all agents are guaranteed not to 
have to spend more battery power than Pmax and they are 
required to contribute to the network connectivity with at 
least the power Pmin- 

The idea of Ref. [7], is to find a trade-off between 
battery consumption and quick and efficient communi¬ 
cation (short network distances) by mixing cheap short- 
range connections and more expensive long-range connec¬ 
tions. The inspiration comes from the small-world network 
model of Ref. [T^ where the authors show that only a 
small fraction of the links need to be long-range—or more 
accurately “any-range”, since they can by construction, 
also can be short—for the entire network to have loga¬ 
rithmically scaling distances (i.e., being efficient, or “small 
world”). In the HKF model, the fraction of any-range links 
is controlled by the parameter q. The construction starts 
from all agents having no links and proceeds as follows. 

1. Pick a random node i with Pi < Pmin- 

2. With a probability q, pick a random node j. If j is dis¬ 
tinct from i and not a member of P^ and if connecting 
i and j would not make Pi or Pj exceed Pmax, then 
add a link between i and j. 
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Fig. 1. An example run of our random walk model of human mobility. In this example we have TV = 10 agents 
and updating time Djj = 30 and the ratio of any-range connections g = 0.1. In the upper panels, we show snapshot 
configurations at t = 20 (a), 28 (b) and 30 (c). These are from the run whose power consumption is displayed at 
the bottom (the times of the snapshots are indicated by triangles). Dashed lines represent range connections that are 
unrestricted in length (at the time of the network construction) and the black line show connections added to nearby 
nodes. The recent trajectories of the nodes are indicated by tails (narrowing with time). 




Fig. 2. An example run of the Levy flight model. The simulation setup and details of the plots are the same as in 

Fig.[g 
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Fig. 3. An example run of the commuter model. The simulation setup and details of the plots are the same as in 
Figs.[^and[^ The periodicity of the commuting behavior 2Dc is the same as the periodicity of the updating Du. The 
commuting behavior is also synchronized so that agents start moving inwards at the time of the updating. 


3. Otherwise (on average a fraction 1 — g of the times), 
attach i to the closest other agent that is not already a 
member of (This could bring the power consump¬ 
tion of i or j above Pmaxj but for a dense enough 
system, probably not much, so we allow such Pmax 
breaches.) 

4. Continue until all nodes have a power consumption 
above Pmin- 

In our simulations, as the agents move, their power con¬ 
sumption can increase over Pmax- We also add a per-agent 
energy cost Eu for the updating itself. The average power 
for a node over one updating period of Du time steps (we 
take time step, without loss of generality, to be one unit 
of time) is thus 


2.3 Mobility models 

Now we turn to our three mobility models. The first one 
is a random walk (RW) with uniform step lengths. This is 
conceptually the simplest possible model and a good null 
model to compare the other models against. We let each 
agent take a step of length 0.05 in a random direction. 
For all three mobility models, we use reflective boundary 
conditions and initialize the agents to uniformly random 
coordinates in the unit square. 

From this simplest mobility model, we generalize in 
two different directions. Our first generalization is as a 
Levy flight (LF), where we use a more realistic distribution 
of the step lengths Z\, namely sampled from a power-law 
distribution [14]. 


P. = 


Du 


Dr. 


P(7 + ^ Pj (t) 


( 2 ) 


The total energy consumption of a node throughout an 
updating period (of duration Du) is the sum of the power 
consumption for all the time steps according to Eq. [l]plus 
Eu- In our simulations we use Pmin = 1 and Pmax = 
2 (just like in Ref. [7]), and Eu = 2. We let P denote 
the power consumption averaged over all agents. We run 
the simulations for 200Du time steps and start taking 
averages after one updating period (to eliminate effects 
of the initial condition). We repeat this procedure n = 50 
times for averages. Unless otherwise stated, we use iV = 50 
agents (or nodes in the networks) and Du = 5. 


P(Z\) cx ® for A e [0.0309...,!]. (3) 

We take the value of the exponent from Ref. [T] (even 
though they study movements of a larger scale than the 
radio links we have in mind). To facilitate comparison, 
we let the LF model have the same average step length 
(0.05), like the other models. We control this by setting 
the lower cut-off to 0.0309... (optionally one could have 
tuned the exponent). The larger cut-off comes from the 
size of the area—we ignore steps that would take us out of 
it. Except the step lengths, this model is identical to RW. 
Second, we mimic the daily routines typical for human 
life mm- We call it the commuter model (CM). We as¬ 
sume agents move with the same step length (0.05) as RW 
but alternately (straight) in toward the center x = y = 0.5 
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or (straight) out toward randomly selected “home” loca¬ 
tions. New home locations are selected at the beginning 
of the simulation (so our results are averages of n sets of 
home locations). They move outwards and inwards over 
intervals of the duration (equaling half the period time) 
Dc- When they are closer than 0.1 from the target (ei¬ 
ther the center or their home location) they take random 
steps of length 0.05 within the circle of radius 0.1 from 
the target. At the start of the simulation (i.e. at the first 
construction of the network), we let the agents be at the 
home location. This will make the effects of Dc more vis¬ 
ible than if one would average the phase difference of the 
two periodic events, or start at some other relative phase. 



3 Results 

3.1 Example runs 

In Fig. [2 we show an example run for the RW mobility 
model. We use N = 10, Du = 30, and q = 0.1. Three 
snapshots of the configurations are shown in panels a, b 
and c. The recent history of the agents is shown by trails 
decaying in width with age. We can see that, by time 
t = 20 and 28 the network has already deteriorated so 
that the power consumption has increased to around two¬ 
fold of what it was after the first construction. You can 
also see the original, short-range links are stretched out 
to the same length as the any-length links. At t = 30, 
the network reconstructed and the power consumption re¬ 
duced from around 2.4 to around 1.4 (if we disregard the 
energy cost Eu of the reconstructing the network, which 
causes the spikes in Fig. &)• After the reconstruction the 
trend of the power consumption is increasing, but there 
are rather large fluctuations due to the small number of 
agents. 

Figure shows the corresponding figure to Fig. for 
the LF mobility model. We can see that, even though the 
fluctuations are larger and more intermittent in the LF 
model, the general evolution of the system is similar to 
the RW model. The behavior for the CM model (seen in 
Fig.§ is, on the other hand, rather different. After an 
updating, the power consumption continues to decrease. 
The reason for this is that, in this case, the commuting be¬ 
havior is synchronized with the updating so that the sys¬ 
tem is updated when the agents are close to their homes. 
After the updating, the agents start moving towards the 
center, which decreases their distances and reduces the 
power consumption. Half the time between the network 
updates, the agents start moving to their home locations. 
Therefore, during this second half of the interval, the dis¬ 
tances and the power consumption increases with time. 


3.2 Effects of the updating time 

Next we turn to the effects of the updating time Du- In 
Fig.i we plot the average power consumption for our 
three mobility models as functions of Du- For the LF and 


Fig. 4. Average power consumption for different duration 
of the updating intervals for the three human mobility 
models. Error bars represent the standard error. The line 
is a fit to the form given by Eq.[^ P « (4.4 ±0.1) ± (9.7± 
0-‘2)lDu- 


RW models, there is an early minimum and a slow in¬ 
crease. The minimum represents the time when the energy 
gain from the optimized configuration balances the cost of 
the updating. For too short Du, one spends too much en¬ 
ergy on constructing the network; for too long Du the 
effects of the gain from the optimized configuration gets 
diluted by the average power-consumption (of what effec¬ 
tively is a random network of the same number of links as 
the one obtained from the HKF model). 

Assume that we can divide an updating interval into 
a time tq when there is an still an power-lowering effect 
from the optimization, so the rest of the time there is no 
other effect of the hkf model optimization than the num¬ 
ber of links. Let Eq be the energy consumed during the 
power-lowering time. Furthermore, let Pu. be the average 
power consumption when there is no correlation between 
the links and the position of the agents. Then, for up¬ 
dating times longer than the optimum, the average power 
consumption is 

p ^ Eu + Eq + Pr{Du - tq) 

Du 

this means P will approach Pr like 


P 


Pr- 


const. 

Du 


(5) 


In practice, we can see this functional for hts reasonable 
for Du > 30 (Fig. |^. For RW one could only achieve a 
reasonable fit with our data for Du > 80 (not shown). 
Once again, this analysis is sketchy (like the assumption 
Tq is a well-defined time scale, the same for all runs). 

The CM model shows a radically different behavior. 
The reason is of course that the power needed to keep 
the radio links operating varies with the position of the 
agents—whether they are closer the center or their home 
locations. Since we (in Fig.|^ construct the network when 
the agents are at the maximal separation, this leads to the 
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sparsest possible network (the agents prioritize keeping 
the connectivity at the expense of performance according 
to the HKF protocol). When the agents move closer to the 
center the network will consume less power, on average, 
than the planned lower limit Pmin- This pattern repeats 
for Djj being a multiple of the periodicity of the movement 
{Dq = 30 in Fig. [^. For higher multiplicities {Du = 60 
and 90) the valleys become broader. We do not have a 
simple explanation for neither this phenomenon, nor the 
minor peaks at other update intervals. If the phase of the 
commuting movement is chosen at random this structured 
pattern gets smoothed out to a curve similar to RW and 
LF. 


3.3 Network structure and optimal updating time 

Next we turn to investigating how the optimal updating 
time depends on the parameters of the HKF protocol. In 
particular, we investigate the dependence on the fraction 
of any-range links q. See Fig.[^ In panel A, we show curves 
(of the same type as Fig. for the RW mobility model 
for some select g-values. For low g, which is the natural 
conhguration, there is a pronounced minimum (as we have 
also seen in Fig. ffl. For larger g-values this minimum be¬ 
comes wider, and eventually any long enough updating 
interval would minimize the average power consumption. 
The reason for the larger g is that the network will be 
sparser (since the HKF protocol prioritizes connectivity 
over robustness and efficient communication). In the ex¬ 
treme situation of g = 1, the locations of the agents do 
not matter at all, so updating the network does not matter 
either. 

Note that this does not mean a large g is necessarily 
preferable. In the present work, we take the HKF model 
as the starting point, which is a model optimizing g for ro¬ 
bustness and performance given an energy budget. These 
quantities are optimized for and intermediate g-values [7] 
and are not changed as agents move after the update. In 
a situation when energy consumption needs to be opti¬ 
mized, and only a minimum robustness and performance 
is required, one would need to optimize g and Djj simul¬ 
taneously. 

In panels b, c and d of Fig. we plot the optimal up¬ 
dating time Du as a function of g. We show the points as 
long as we can locate the minimum with a precision higher 
than five time units. As the minimum gets flatter the im¬ 
precision of locating it explodes. To highlight this, we plot 
the standard deviation, rather than the standard error, in 
panels b-d. We show the region of ill-defined minima— 
where the valley is so flat compared to the fluctuations, 
so that we cannot pinpoint the minimum—as shaded in 
panels b, c and d. For the RW and LF models, the optimal 
updating intervals grow with g. This effect is stronger for 
RW than LF (which is in contrast to the Uiy-dependence 
of P, Fig. 1^ where LF is most sensitive to parameters 
than RW). The optimal updating interval of CM model 
increases very slowly. It is fair to call it independent of g. 
We do not have any simple mechanistic explanation for 




q q q 

Fig. 5. The updating duration minimizing the power con¬ 
sumption as a function of the fraction of any-range links. 
In the shaded areas there is no intermediate optimum— 
the cost of updating the network outweighs the benefits 
(so longer Du is better). The error bars of panel a rep¬ 
resent standard errors, while in b-d they represent the 
standard deviation. 



Fig. 6. Average optimal power consumption as a function 
of the fraction of long-distance connections g for the three 
models of human mobility. Error bars represent standard 
errors. 

this, but probably Du is pinned to the periodicity of the 
commuter movement. 


3.4 Network structure and optimal power consumption 

Our final analysis concerns the average power consump¬ 
tion at the optimized updating intervals as a function of 
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the fraction of any-range links. The results are shown 
in Fig. These curves have maxima at intermediate q- 
values. A reason for the increase at small q-values is that 
during the network construction process, one is allowed to 
use more short-range links which means one can reach over 
Tmin without overshooting much. Allowing long-range links 
means that the lastly added link of several nodes is long- 
range which could bring the power consumption close to 
Tmax and thus increase the average optimal power con¬ 
sumption. For large ^-values the effect, as mentioned, is 
that the network becomes sparser which—as soon as the 
memory of the original configuration has deteriorated by 
the movement of the agents—lowers the average power 
consumption. At q = 1, the only difference between the 
LF and RW is that the larger fluctuations of LF deteri¬ 
orates the optimization faster (so LF’s value becomes a 
bit higher). The power consumption for the CM model is 
lower than both RW and LF since, it the links shrink when 
the agents are close to the center and the DC is a multiple 
of Du, so that the network construction takes place when 
the network is as spread out as possible. This effect would, 
of course, change if we initially constructed the network 
at some other phase of the agents’ movement. 


4 Discussion 

In this paper, we have explored the impact of the move¬ 
ment patterns of agents on the power consumption of wire¬ 
less ad hoc networks. This is a question about engineering 
emergent phenomena and as such it is connected to many 
fields—electric engineering, the theory of distributed com¬ 
puting, and aspects of the physics of complex systems such 
as small-world networks and Levy flights. Roughly speak¬ 
ing our objective is to simultaneously achieve five proper¬ 
ties: 

1. Low global power consumption. 

2. A fair distribution of the power consumption. 

3. A short hopping distance. 

4. A network robust to failures and attacks. 

5. Parsimony. 

In a practical implementation, other issues need to be 
added to this list, like the computational efficiency and ro¬ 
bustness of the algorithms M. Of the five points above, 
points 3 and 4 only relate to the static network topol¬ 
ogy, so we focus on I in this work. (While points 2 and 

5 are directly taken care of by the network-construction 
protocol.) 

We find that the movement patterns do affect when the 
network need to be updated in non-trivial ways. The net¬ 
work need to be updated most rarely with the RW mobil¬ 
ity model, and most frequently with the LF model (Fig.[^. 
In other words, heterogeneities of the step lengths induce 
a need for frequent updating (even though the agents have 
the same average speed in all simulations). Both RW and 
LF have increasing optimal updating intervals as func¬ 
tions of the fraction q of any-range links (that tunes the 
balance between robustness and efficient communication 


in the network). The CM model, on the other hand, has 
the same optimal updating interval regardless of q. An¬ 
other finding (Fig. [ 6 ) is that the most efficient power con¬ 
sumption of both the RW and LF mobility models have 
peaks as functions of q. Interestingly, the minimum aver¬ 
age power consumption happens at opposite extremes for 
RW and LF (9 = 0 for RW and q = I for LF). Also for 
this analysis, the CM curve is flatter than the others. 

There are, of course, many other directions to explore 
in such a complex problem as the interaction between ad 
hoc wireless network and mobile agents. Perhaps the least 
realistic assumption in our work is that any pair of devices 
can establish a radio link throughout the studied area. In 
a situation where the agents could move outside of the 
range of their devices, the overall power consumption (as 
in Fig.[^ would not be the most critical performance met¬ 
ric. One would also need to evaluate the connectivity of the 
network, or the average graph distances (cf. Ref. [H]). For 
RW and LF, these functions are monotonically decreasing 
and the overall picture of our analysis would be the same. 
For CM there could be a situation where the connectivity 
and distance metrics oscillates with the periodicity of the 
movement. We could add another layer of realism by al¬ 
lowing the population to fluctuate. This would, however, 
probably not alter our conclusions, at least if the connec¬ 
tions are not so expensive that the connectivity of the 
entire system is an issue. 

We think there are many interesting future questions 
relating to the present work. Testing other mobility mod¬ 
els would be one promising direction. One could for ex¬ 
ample assume the agents’ movement was affected by their 
role in the network. If agents tried to (selfishly) mini¬ 
mize their power consumption, or (unselfishly) optimize 
the communication efficiency and robustness of the net¬ 
work, that could give rise to a yet more complex behavior 
than we have seen. Another interesting extension of our 
work would be to use empirical data for human mobil¬ 
ity mm instead of models of the mobility, or adapt it 
for vehicular traffic [mg. 
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